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Abstract 

We examine the contributions of odd-parity nuclear operators to the two-neutrino double 
beta decay 0+ 0+ amplitude, which come from the P-wave Coulomb corrections to the 
^ ! electron wave functions and the recoil corrections to the nuclear currents. Although they are 
'sj" . formally of higher order in or f /c of the nucleon than the usual Fermi and Gamow- 

-p^ ■ Teller matrix elements, explicit calculations performed within the QRPA show that they are 
significant when confronted with the experimental data. 
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The two neutrino double beta decay {2v(3(3) 0^ 0^ is the rarest process observed 
so far in nature. As such, it offers a unique opportunity for testing the nuclear physics 
techniques for half-lives Tfj^ ~ (10^^ - 10^3) years % |, |, |, |, ||. But, before resorting 
to a particular nuclear model, the expression for the 2v(3(3 decay rate has to be figured 
out from the weak Hamiltonian. Except for the work of Williams and Haxton [H, the 
allowed approximation (AA) has been universally adopted in all previous theoretical studies. 
This approximation is obtained in the long-wavelenth limit for the outgoing electrons and 
neutrinos by retaining only the Fermi (F) and Gamow- Teller (GT) operators in the non- 
relativistic impulse approximation for the hadronic current 

J'ij) = E [i9v - 9A^n ■ ^n)g'' + (^7^^^ " gvvt)g''] 5{r - r^. (1) 

n 

Here the operator r+ transforms a neutron into a proton and cr„, r„ and v„ stand, respec- 
tively, for the spin, the position and the velocity of the nucleon n. The well known result for 
the half-live is ||, |Tg, [TT], 
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Tf^Q+ - 0+)J = ^2.(W^o,^)|-M2.r, (2) 
where QiWo, Z) is the phase space kinematical factor, and 

(0^k+|0+)(0+|r+|0+) 



■M2U = J^2U - -9V 2^ ^ ^ , iTT/ 



is the nuclear matrix element in the AA, with all the notation having the usual meaning 



10, 11 



The retardation of the 0+ 0+ 2vj3(3 transition rates mostly comes from the phase 

space factor. Yet, the nuclear correlations also slow down these processes. Their effect can 

be estimated by comparing the measured |A^2i/| values, shown in table with the naive 

^The 2i//3/3-decay is a second order nuclear weak transition, analogous to the electromagnetic nuclear 
double gamma decay. It could be interesting to notice that the half-life for the 77 0^ Q'l transition in 
9°Zr is only tIJ^ - 1.4 x IQ-"^ seconds §. 
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estimate, Ai2i/ ~ —0.3 that is obtained within the closure approximation for an effective 
axial charge g'^-^ = f 0. The resulting suppression factor ranges from 1/2 for ^°°Mo to 1/15 
for i30Te. 

The quasiparticle random phase approximation (QRPA) is the nuclear structure method 
most widely used to deal with the quenching of the moment Ai2i'- We know that in this 
model the allowed moment Ai2u is very sensitive to the ground state correlation (GSC) 
within the particle-particle (PP) channel. In a recent work it has been shown both that: 
(a) the restoration of the isospin symmetry leads to A^^^(J'' = 0"*") = 0, and (b) the exact 
QRPA calculations for 7W^^(J'^ = l"*"), when evaluated with a zero-range force, can be nicely 
fitted by the formula 

Mi{r = 1+) = Mi{r = 1+; t = 0)[^, (4) 

where t = ff^/f^"*'^ is the ratio between the spin-triplet strength in the PP and the spin- 
singlet strength in the pairing channel. The free parameters to and ti denote the zero and the 
pole parameters to and ti, with ti ^ to, denote the zero and the pole of A^^^, respectively. 
In the same work it has been suggested that this result is of general validity and that 
any modification of the nuclear hamiltonian or of the configuration space cannot lead to 
a different functional dependence. Several alterations of the QRPA have been proposed 
that might change that behavior. They include higher order RPA corrections |M, nuclear 



deformation |T^, single-particle self-energy BCS terms |]T6[ and particle number projection 
r^ . Yet, none of these amendments inhibits the matrix element pass through zero 



near the "natural" value of the PP strength {t ~ 1.5). f\ 

We also know that the GSC do not affect the Oi//5/3-decay probability in the same way. 
The reason for that is that the exchange of the virtual neutrino in the Ou/SP process gives 
rise to the neutrino potentials, which depend on the distance ri2 = |ri — between the 
two vertices where the single /3-decays take place. If one expands this dependence on ri2 in 



^When a physical quantity has a zero (or near zero) a conservative law should, very likely, be its origin. 
Thus, the behavior of A^^^, described by has been attributed to the partial restoration of the SU(4) 
Wigner symmetries This is not surprising since the operators r^cr are infinitesimal generators of the 
SU{4:) group. 
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multipoles depending on ri and r2, one sees that, besides exciting F and GT states, one is 
also going through virtual states with spin and parity J'^ 7^ 0"*", 1^. The GSC are not so 
important for the forbidden transitions {L > 0) and thus they are not so drastically reducing 
the decay rate for the Oz//3/5-decay as for the 2z//?/3-decay. In fact, it turns out that for t ~ 1.5 
the dominant contribution to Aiou comes from multipoles with L > ||T3 |. 



Motivated by the above arguments and by the hope that some favorable nuclear structure 
physics could enhance the higher order contributions to A^2i/, we study here the effects of the 
2vl3(3 matrix elements with L = 1. It will be assumed that the Coulomb energy of the electron 
at the nuclear radius is larger that its total energy. This leads to the ^-approximation, which 
has been extensively utilized in the study of the first forbidden (FF) single 13 transitions ||19|| . 
Besides, it is well known that, when this approximation is applied, the FF observables can 
be expressed exactly as the allowed ones except for the replacements p 

r+gy _ -r+X; t^qact ^ -r+Y, (5) 

with 

Y = gvv + ^{gvir - gAcr X r), (6) 

where (in natural units) 

e = ^ = 1-18 ZA-'/\ (7) 

Z and R being the charge and the radius of the daughter nucleus. Q This result is obtained 
by retaining the velocity dependent terms in the hadronic current (|I]) and by considering 
only: i) the 5*— wave for the neutrino wave function 4''^{s), and ii) the S— wave and the C, 
correction term in the P— wave for the electron wave function il''^{r, s), i.e.. 




■^We discuss here only the /3 -decay. For positron emission the following changes are made: Z — > —Z and 
9 A -gA- 
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and 

s) = ^ ^"(^■^)(^ + '") 11 H. ,(f^ . r)| ( ^ ) (9) 

where Fo{Z,e) is the Fermi factor and all the remaining notation is that of ref. 

Before proceeding let us recall that the contributions of all five FF moments are of the 
same order of magnitude. First, in the extreme single particle model 

= ^fi^f + 1) - ^/(^/ + 1) - ^'^(^'^ + 1) + ^^(^^ + (10) 

Next, the velocity dependent /3-moments are related in a simple way to the corresponding 
position operators, when the harmonic oscillator radial wave functions are employed. One 



gets 19, 20 



T7V- TV = \-\ = ^/ ~ 11 

where N and ujq = 80 A~^^'^ stand for the principal quantum number and the oscillator 
frequency, respectively. The most relevant single-particle transitions, in any nuclear model 
calculation, are those with Ni = Nf + 1 (see e.g., table 2 in ref. 12^). Consequently the 
contributions of the operators cr ■ v and v always add destructively to those coming from the 
operators of iC,(T ■ r and i^r (see eq. (^). In addition, the ratio of the corresponding matrix 
elements is ~ —68/Z, i.e., of order of unity for medium heavy nuclei. 

The conserved-vector-current (CVC) constraint leads to an alternative expression for the 
total matrix element < v > |^ 

< V >= -(2.4^ + Eo) < ir >, (12) 

where Eq is the energy difference between the initial and final states. (In principle, this 
relation is exact as it includes exchange effects, induced interactions, etc. [|l9l-) Besides, it 
is now well established [2^ that the operator cr • v is enhanced by the one-pion exchange 



current contributions in nuclear medium pi]. Inspired by this effect, which can be accounted 



''Note that the operators icr ■ r (ir) and cr • v (v, cr x r) have different Hermitian conjugation properties 
[expressed in equations ([Io| ) and (^l])]. This makes that in a RPA calculation the signs of the ratios ( |l0[ ) and 
are different for the forward and backward going contributions. 
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for by renormalizing the corresponding matrix element by the factor 1.7 g^'^'^ = 2.0 (see also 
ref. ||21[]), Williams and Haxton |^ have pointed out that the moment < <t ■ v > could 
have a significant effect on the rates for 2z//3/3-decay. They also performed calculations in 
the Nilsson pairing model with the quadrupole-quadrupole core-polarization effect included. 
The result was to increase the AA predictions for of ^^Ge, ^"^Se, ^^^Te and ^^^Te nuclei 
by 15%, 25%, 15% and 20%, respectively (see also ref. |^). Yet from the above discussion 
it is clear that the consequences of the remaining four FF moments have also be examined 
simultaneously. | Such an examination is the main topic of this paper. 

The expression for the half-life for the 0^ — O"*" 2z//?/5-decay, within the ^-approximation, 
has been derived following the procedure of Doi et al. and using the relations (H) and 



(P) for the lepton wave functions. A tedious but straightforward calculation ||26[ leads to the 
equation for the half-life that can again be cast in the form (P, with the same kinematical 
factor and the transition amplitude 

M2u^Ml = Mi + M^f, (13) 

where 

j^FF^ _ j^ {0t\r'-X\0-^){0-^\T+X\0t) 
(0+|r+Y|l^)-(l]v|r+Y|0+) 

is the contribution coming from the FF moments. It is also important to note that the shape 
of the electron spectrum is the same as in the AA. 

As the spin singlet components predominate in the ground state of the even-even nuclei, 
the two terms in (^ and ( p^lf ) have a tendency to sum up coherently. On the other hand, 
it is self-evident from relations (|]) that, independently of the nuclear model employed in 
the evaluation of i ^^e role played by the virtual states with J'^ = 0~ and J'^ = 1^ 
critically depend on the way in which the individual matrix elements combine to build up 



^Williams and Haxton g have also hinted that, besides < cr-v >, other parity-forbidden matrix elements 
should be considered. Particularly, because of the CVC relation, they have pointed towards the importance 
of < V >. 
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the total matrix elements of X and Y. Coherent combinations would strongly enhance the 
effect of M.2if- Yet, from the above discussion we have learned that the moments < cr ■ v > 
and < icr ■ r >, as well as < v > and < ir >, always tend to add destructively. Hence the 
magnitude of M.2if hangs crucially on the matrix element < cr x r > (both on its magnitude 
and its sign relative to the sign of < v >). 

The numerical results for the matrix elements M.2u presented in table 0. They 
were obtained through the procedure adopted in our previous works |]18|, i.e., by using 
the same residual interaction, the same configuration space, etc. The moments labeled 
A4f/(pair) were evaluated in the independent quasiparticle approximation, and those indi- 
cated by M.2v{Q^P^) within the QRPA. The matrix element < cr- v > has been estimated 
by using the relation (pi]), properly renormalized by the one-pion exchange currents. The 
CVC relation ([T^) has been used in the evaluation of the moment < v >. The particle- hole 
interaction reduces the forbidden moments by roughly a factor of 2. But, at variance with 
the moments M.2v^ which are widely uncertain in the QRPA (due to the proximity of ti to 
to in eq. (^), the effect of the particle-particle interaction on M.2v is always very small. For 
all the nuclei, the contributions of the virtual states with J'^ = 0~, shown parenthetically 
in table 0, are smaller than those coming from the states with J'" = 1~ . There are two 
reasons for that. First, the individual single-particle matrix elements of the operator Y are 
usually larger than those of the operator X. This is caused mainly by the destructive inter- 
ference between the individual FF matrix elements mentioned above. Among the former the 
most important are those with (/|cr x r\i) = {f\ir\i), in which case {f\Y\i) = (/|v|i) (for 
g'^^ = 1). Second, for both type of transitions the partial contributions show a pronounced 
coherence, but the configuration space for the 1~ is bigger. For the sake of comparison, the 
matrix elements A^^j, are displayed in table 2 also. 

We note that the calculation done by Williams and Haxton IQ yields considerably larger 

^It is worth noting that A^f^T can be derived from by performing the substitution (||) in i.e., in 
the analogous way that are obtained the FF observables from the aUowed ones in the simple /? decay. This 
common future of the single and double beta decays came somewhat as a surprise. One should remember 
here that, because of the antisymmetrization in the momentum and spin-orientation, the 2v(3l3 amplitude 
is not just a simple product of two single (3 decays. Yet, a careful examination of the analytic calculations 
shows that it is precisely the anti-symmetrization procedure that establishes such a correspondence. 
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contributions of the = states that our QRPA evaluation. For instance, for the 
s^Se ^^Kr and ^^°Te ^^'^Xe transitions they get M^fiO') ^ 0.010. The dis- 
crepancy with our results can be attributed only partially to the effect of the matrix element 
< icr ■ r >, which has not been considered in ref. P]. By neglecting the last one we obtain 
that M^f{0-) = 0.006 for both nuclei. The remaining discrepancy has to be attached to 
the difference in the nuclear models employed in ref. and here. 

From the results shown in table |] it is clear that in the QRPA the matrix elements 
A^^jf are significant when compared with the A42u moments. But, the real relevance of 
■^2if comes from the confrontation with the experimental data. In fact, we see that they 
are roughly comparable to \ A42u\{Gxp) needed to explain the half lives of ^^^Te, ^^°Te, and 
^^^Nd nuclei. For the remaining four nuclei displayed in tables |l| and ^ the moments A^^jf 
are relatively small in comparison with the measured amplitudes. Yet, even in this case the 
parity forbidden moments might become significant in the theoretical evaluation of the 2v(3l3- 
decay half-lives, as they have a tendency to cancel against the allowed moments. Besides, 
the sensitivity of 2vf3l3 experiments has been greatly improved in recent years, and there is 
no reason to believe that further improvements are not still possible. 

In closing, we point out that there are other forbidden terms that should be studied. 
For instance, we are in process of examining the contributions arising from virtual states 
J'^ = 2~, which might modify both the electron spectrum shape and the Q{Wo,Z) factor 
predicted by the AA. This could be of relevance for the excess of high-energy electrons seen 
in the (3(3 spectra and presently ascribed to the majoron emission. Finally, we notice that in 
the new class of majoron models the /3/3-decay proceeds via the odd-parity nuclear operators 

H • 

We warmly thank Dubravko Tadic for fruitful discussions. 
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Tables 



Table 1: The measured half-lives ^"^^ phase space function G2u{W{)-, Z), which were 
used to extract the values of |A^2i/|(exp), given in natural units. 



Nucleus 






-1 


\M2u\iexp) 




14.3 ±0.17'^^ 


5.39 10" 


-20 


0.114 ±0.007 




1 nS+O-26 b) 

l-Uo-o.oe 


1.80 10" 


-18 


n7i+°-°°2 

U-U 1 i-0.007 




0.39 ± 0.09 


7.66 10" 


-18 


0.058tnS^ 




n 1 1 c;+0.030 d) 
U.iiO_o.o20 


3.91 10- 


-18 


0.1491°:°!^ 




(7.7 ±0.4) 10^") 


3.53 10- 


-22 


0.019 ±0.001 




27 ± 1 


1.98 10- 


-18 


0.014 ±0.001 




n 1 7+0.14 /) 
U--L'-0.09 


4.91 10- 


-17 


U.UOd_o.o09 



") (laboratory data) ref. JIl 
^) (laboratory data) ref. 
'^) (geochemical data) ref. 
'^) (laboratory data) ref. [Q 
*^) (geochemical data) ref. 0] 
■^) (laboratory data) ref. ||^ 
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Tabic 2: Calculated allowed and first forbidden matrix elements (in natural units) for an 
effective axial charge g'^-^ = 1. The column labeled Ai2y {pair) is evaluated using the inde- 
pendent quasiparticle approximation, while the column A1f/(Qi?PA) includes the effects of 
both the particle-hole and the particle-particle interactions. The values in parenthesis stand 
for the contributions of virtual states with J'^ — 0~ to the total moments. 



Nucleus 




(pair) 




iQRPA) 


MliQRPA) 




-0.015 


(-0.005) 


-0.008 


(-0.004) 


0.050 




-0.019 


(-0.006) 


-0.009 


(-0.004) 


0.060 




-0.033 


(-0.010) 


-0.014 


(-0.006) 


0.010 




-0.035 


(-0.010) 


-0.014 


(-0.006) 


0.051 


128yg 


-0.020 


(-0.004) 


-0.012 


(-0.003) 


0.059 




-0.019 


(-0.004) 


-0.012 


(-0.002) 


0.048 


150^^ 


-0.060 


(-0.009) 


-0.031 


(-0.008) 


0.033 
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